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A method for calculating the density of defects to be expected from Compton 
electrons formed by incident y-irradiation of a given material is described. The 
calculation involves experimental knowledge of the rate of introduction of defects, 
d’V d /d0, where N d is the density of defects and 4 > the integrated electron flux, the 
range-energy curve for electrons, and the theoretical spectrum of Compton 
electrons given by the Klein-Nishina formula. A specific calculation has been 
made for the density of defects at the band energy level 0.17 eV below the conduc- 
tion band energy, E c , inn-type silicon irradiated by 1.37 MeV Co 60 y- rays. 

Sonder and Templeton ( found a defect production rate of 1.35 x 10” 3 defects per 
cm path of Co 60 photons in 2 Ohm-cm silicon containing oxygen. Using this 
production rate, they calculated the defect production rate per Compton electron, 
which turned out to be by a factor of ten smaller than the experimental produc- 
tion rate for accelerator electrons of comparable energies. The result of the 
present calculation is in reasonable agreement with their experimental data of 
1.35 x 10 -3 cm’ 1 . It corrects their theoretical calculation by determining the 
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equilibrium spectrum of the Compton electron flux, N(E), which results from the 
degradation of electron energy in ionizing collisions. 

The first step is to calculate the original Compton spectrum^ 2) as cut-off 

J 

by the maximum electron energy, which is proportional to d<t> c /de where C> c is 
the Compton cross section per target electron which acquires an initial energy e. 
This spectrum is shown in Figure 1. 

It is then possible to calculate the number of equilibrium electrons, N(E), 
with energy between E and E + dE crossing a unit area within an infinite medium 
by summing over all electrons which will have the energy E at the plane located 
at z = 0. For a unit fluence of gammas in the z-direction, 


r E max /.O) /•n/2 

N(E)dE = de dz 
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where <t> is the azimuthal angle for the z-direction, N e is the number of electrons/ 
cm 3 and P(z, 9, e, E) dE is the probability that an electron of energy e which 
has been scattered at an angle 9 with respect to the incident direction of the 
gamma ray at a distance z away from the flux plane will have an energy between 
E and E + dE. P(z, 6, e, E) will on the average be unaffected by elastic scatter- 
ing of the electrons by nuclei since as many electrons are scattered away from 
as toward a given polar direction 9. The expression for P, n'ormalized over all 
energies, is 


P(z, 9, e, E)dE = -gf- | E 1 S |R(e) - R(-E) -.z/cos 9 1 dE 
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where S is the Dirac delta function, R is the range, and E’ is given by: i 

R(E’ ) =R(e) - z/cos 8. The expression for d<£> c /de(e, 8, <£), normalized over i 
all angles, is: 

t 

: d<D c /deS (8 -9\ 

d<l> c /de (e, 6, 4>) = “ 2~ s ~ in -q € ' (3) 

where 0 e is the angle of scattering corresponding to e. The indicated integration 
gives: • 

fjn f E -.x dO. 

N(E)dE = dE N e -gg- | E de — cos 6 t . (4) 
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The final integration was performed numerically and the range-energy relation 
for silicon was assumed equal to the expression for aluminum when the density, 
P, for silicon was used. Thus, from Marton (2) , 

R(cm) = (0.412/,o)E ( 1 - 265 ' 0 - 09S41nE ) 

when E is given in MeV and p in gm/cm 3 . The equilibrium spectrum is also 
given in Figure 1. j 

A graph of dN/dd> versus E is shown in Figure 2 for the defect located at 
E c - 0.17 eV. The curve represents the average of the latest data by Carter and 
Downing (3) , which appears to be the most extensive and consistent data available. 
The data by Flicker, Loferski, and Scott-Monck (4) , representing relative values 
only is normalized to the Carter and Downing data at 800 KeV. It is realized 
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that the curve may have a spread of a factor of two if only because of variations 
in dN/d<f> which occur for samples of different resistivities. 

The final calculation for N d is given by 

N d = dE N(E) (dN d /d<£) (5) 
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or N d = 2.2 x 10 _3 /cm 3 per unit integrated flux of y's. This number is a factor 
of two larger than the experimental data of Sonder and Templeton ( . However, 
the latter authors have estimated that their value should be low by 30-50% due to 
the geometry of the counting device. 
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